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Abstract 

The  interaction  of  atoms  and  molecule  with  surfaces  is  of  great  technological  relevance.  Both 
advantageous  and  harmful  processes  can  occur  at  surfaces.  If  an  atom  or  molecule  impinges 
on  a  surface,  it  can  either  scatter  back  into  the  gas  phase  or  become  adsorbed  on  the  surface. 
Molecules  can  furthermore  undergo  chemical  reactions  at  the  surface.  All  these  processes  are 
accompanied  by  energy  transfer  between  the  impinging  projectile  and  the  substrate. 

The  simulation  of  the  dynamics  of  the  gas-surface  interaction  still  represents  a  considerable 
challenge  since  the  coupling  of  a  low-dimensional  object,  the  impinging  atom  or  molecule,  to  the 
substrate  with  in  principle  infinitely  many  degrees  of  freedom  has  to  be  modeled.  Furthermore, 
depending  on  the  mass  of  the  atom  or  molecule,  quantum  effects  both  in  the  molecular  motion 
as  well  as  in  the  excitation  of  the  substrate  have  to  be  taken  into  account.  In  this  lecture,  the 
quantum  and  classical  methods  required  for  the  simulation  of  gas-surface  dynamical  interactions 
will  be  reviewed.  Furthermore,  the  main  processes  occuring  in  the  interaction  of  atoms  and 
molecules  with  substrates  will  be  illustrated  using  quantum  calculations  and  classical  molecular 
dynamics  simulations. 


1  Introduction 


Understanding  the  interaction  of  atoms  and  molecules  with  surfaces  plays  an  important  role  in  a 
wide  range  of  technologically  relevant  applications  [1].  Among  those  are  the  heterogenous  catalysis 
-  the  majority  of  reactions  in  the  chemical  industry  employ  catalysts;  crystal  growth,  which 
determines,  e.g.,  the  quality  of  semiconductor  devices;  corrosion  and  lubrication,  which  influences 
the  durability  of  mechanical  systems;  or  friction,  which  determines  the  energy  transfer  to  surfaces, 
as  for  example  in  the  interaction  between  the  atmosphere  and  rockets  or  high-speed  missiles  which 
causes  significant  temperature  rises  that  lead  to  chemical  reactions  and  particle  ionizations. 

In  order  to  study  the  interaction  of  molecules  with  surfaces,  it  is  mandatory  to  have  a  reliable  de¬ 
scription  of  the  interaction  energetics.  The  basic  quantity  reflecting  this  interaction  is  the  potential 
energy  surface  (PES)  that  corresponds  to  a  hyperplane  in  the  multidimensional  configuration  space. 


GroB,  A.  (2007)  Simulation  of  Gas-Surface  Dynamical  Interactions.  In  Experiment,  Modeling  and  Simulation  of  Gas -Surface  Interactions 
for  Reactive  Flows  in  Hypersonic  Flights  (pp.  4-1  -  4-26).  Educational  Notes  RTO-EN-AVT-142,  Paper  4.  Neuilly-sur-Seine,  France: 
RTO.  Available  from:  http://www.rto.nato.int/abstracts.asp. 
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Figure  1:  Schematic  illustration  of  the  role  of  a  catalyst  using  a  two-dimensional  representation 
of  the  potential  energy  surface.  A  catalyst  provides  a  detour  in  the  multi-dimensional  PES  with 
lower  activation  barriers 

Nowadays  potential  energy  surfaces  can  be  mapped  out  in  great  detail  by  first-principles  electronic 
structure  calculations,  typically  based  on  density  functional  theory.  From  the  PES,  quantities  such 
as  adsorption  or  binding  energies,  transition  state  configurations  and  energies  or  minimum  energy 
paths  can  be  derived. 

Potential  energy  surface  can  also  be  used  to  illustrate  processes  at  surfaces.  An  example  is  given  in 
Fig.  1  explaining  the  way  a  catalyst  works.  The  presence  of  a  catalysts  lowers  the  activation  barrier 
for  a  particular  chemical  reaction,  however,  this  involves  a  detour  in  the  multi-dimensional  PES  on 
the  path  from  the  reactants  to  the  products.  In  heterogeneous  catalysis,  this  usually  means  that 
the  reactants  have  to  be  adsorbed  on  the  catalyst  surface,  where  the  activation  barrier  is  much 
smaller  than  for  example  in  the  gas  phase.  Hence  the  reaction  rate  is  enormously  enhanced  in  the 
presence  of  a  catalyst  since  the  rate  depends  exponentially  on  the  barrier  height. 

For  a  true  understanding  of  processes  on  and  at  surfaces,  however,  the  static  information  from 
total-energy  calculations  is  often  not  sufficient.  Furthermore,  experiments  do  usually  not  yield 
direct  information  about  the  potential  energy  surface,  but  rather  determine  reaction  probabilities 
or  transfer  rates.  In  order  to  calculate  these  quantities  and  allow  for  a  genuine  comparison  between 
theory  and  experiment,  dynamical  or  kinetic  simulations  are  required.  Unfortunately,  processes 
such  as  chemical  reactions  at  surfaces  often  consist  of  many  elementary  steps  which  are  too  complex 
to  be  studied  as  a  whole.  Therefore  in  surface  science  one  tries  to  understand  reaction  mechanisms 
by  breaking  them  up  into  simpler  steps  which  are  then  studied  under  well-defined  conditions  [2] . 

In  this  chapter,  I  will  briefly  review  the  theoretical  methods  necessary  to  determine  the  dynamics 
of  processes  at  surfaces.  After  presenting  classical  and  quantum  methods,  scattering  at  surfaces 
will  be  addressed.  A  substantial  fraction  of  the  chapter  will  then  be  devoted  to  the  discussion  of 
atomic,  molecular  and  dissociative  adsorption.  Finally,  the  system  02/Pt(lll)  will  be  discussed  as 
an  example,  where  scattering,  molecular  and  dissociative  adsorption  can  all  occur. 
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2  Molecular  Dynamics 


In  order  to  follow  the  temporal  evolution  of  a  dynamical  system,  the  corresponding  quantum  or 
classical  equations  of  motion  have  to  be  solved.  They  can  be  derived  from  the  basic  equations 
describing  the  interaction  of  atoms  and  molecules  with  each  other,  namely  the  non-relativistic 
Schrodinger  equation  on  which  most  of  chemistry  and  solid-state  physics  is  based.  It  is  of  the 
general  form 

tftf({R,r})>  =  £*({R,r})  .  (1) 

where  the  R  are  the  ionic  coordinates  and  r  the  electronic  coordinates.  It  is  well-known  that 
a  complete  analytical  solution  of  the  Schrodinger  equation  taking  into  account  both  ionic  and 
electronic  degrees  of  freedom  is  not  possible  except  for  simple  cases.  One  common  approach  is 
to  assume  that  -  because  of  the  large  mass  difference  between  electrons  and  the  nuclei  -  the 
electrons  follow  the  motion  of  the  nuclei  adiabatically.  This  is  the  so-called  Born-Oppenheimer 
approximation  [3]. 

In  practice,  in  the  Born-Oppenheimer  approximation  the  nuclear  coordinates  are  fixed  so  that  they 
enter  the  Schrodinger  equation  no  longer  as  variable  but  just  as  parameters.  Thus  one  obtains  an 
electronic  Schrodinger  equation 


He  1  ({R})  VKM)  =  Ee  1  ({R})  |<KM). 


(2) 


The  many-electron  ground  state  energy  Fiei({R})  then  defines  the  potential  energy  surface  for  the 
motion  of  the  nuclei.  For  extended  systems  the  most  efficient  approach  to  determine  the  many- 
electron  ground  state  energy  Fiei({R})  from  first  principles  is  density  functional  theory  (DFT)  [4] 
in  combination  with  the  supercell  concept.  Once  the  electronic  ground  state  energy  is  obtained,  it 
can  be  plugged  into  the  Schrodinger  equation  for  the  nuclei, 


Vr,  +  Ee  i({R}) 


En  ucl  $({R}), 


(3) 


where  Enuc\  is  now  the  energy  relevant  for  the  dynamics  of  the  nuclei.  Alternatively,  the  potential 
energy  surface  can  be  used  to  solve  the  classical  equations  of  motion, 


, ,  c>2  Q 

1  m 


Eel{{R-m}) 


(4) 


Instead  of  Newton’s  equation  of  motion,  also  Hamilton’s  equation  of  motion 

dH  .  dH 


(5) 


can  be  solved. 


The  solution  of  the  equations  of  motion  can  be  obtained  by  standard  numerical  integration  schemes 
like  Runge-Kutta,  Bulirsch-Stoer  or  predictor-corrector  methods  (see,  e.g.,  [5]).  Very  often  the 
rather  simple  Verlet  algorithm  [6,  7]  is  used  which  is  easily  derived  from  a  Taylor  expansion  of  the 
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trajectory. 


r  i(t  +  h) 


dYi 

h 2  d?Ti 

h 3  d?Yi 

dt 

h= 0  2  dt 2 

h= 0  1  6  dt3 

r i(t)  +  h 


.  h?  Fj(t) 

r  i(t)  +  hvi(t )  +  — - 

2  m 


+  ... 


h3  d3r i 

+  YW 


h= 0 

+  ... 


h= 0 


(6) 


Here  we  have  introduced  the  velocity  v,;  =  rj.  Furthermore,  we  have  used  Newton’s  equation  of 
motion  to  include  the  force  F,;  acting  on  the  i-th  particle.  Analogously  we  can  derive 


.  , .  .  .  ,  .  .  h?  FAt)  h3  d3Ti 

r,((  -  h)  =  r,(f)  -  h  v,(f)  +  -  —  -  ^  w  + 


h= 0 


Adding  (6)  and  (7)  yields  the  Verlet  algorithm  [6] 

r i(t.  +  h)  =  2r i(t)  -  r i(t  -  h)  +  h2  +  0(h4)  . 


m. 


(7) 

(8) 


The  accuracy  of  the  numerical  integration  of  the  equation  of  motion  can  be  checked  by  testing  the 
energy  conservation.  In  order  to  evaluate  the  kinetic  energies,  the  velocities  at  time  t  are  needed. 
Note  that  they  do  not  explicitly  appear  in  Eq.  (8).  They  can  be  estimated  by 


v,(t)  =  (9) 

However,  the  kinetic  energy  evaluated  with  Eq.  (9)  belongs  to  the  time  step  prior  to  the  one 
used  for  the  positions  (8)  which  enter  the  evaluation  of  the  potential  energy.  This  problem  can  be 
avoided  in  the  so-called  velocity  Verlet  algorithm  [7] 


r  i(t  +  h) 
v*(t  +  h) 


r i(t)  +  h  Vi(t )  + 


Fj(t) 
2  m 


Vj(t)  +  h 


F  i(t  +  h)  +  F  i(t) 
2m 


which  is  mathematically  equivalent  to  the  Verlet  algorithm. 


(10) 


In  Fig.  2,  various  processes  that  can  occur  when  atoms  or  molecules  are  impinging  on  a  surface 
are  illustrated.  For  all  atoms  that  are  heavier  than  hydrogen  or  helium,  the  quantum  effects  in  the 
dynamics  are  often  negligible  [8].  Hence  molecular  dynamics  simulations  involving  the  solution  of 
the  classical  equations  of  motion  are  an  appropriate  tool  to  determine  the  time  evolution  of  most 
chemical  systems.  Even  if  hydrogen  is  contained  in  the  system  as  is  generally  the  case  for  organic 
molecules,  still  the  results  of  classical  dynamical  calculations  might  be  meaningful  qualitatively 
or  semi-quantitatively.  However,  for  a  real  quantitative  description,  a  quantum  treatment  is  nec¬ 
essary.  Furthermore,  there  are  certain  phenomena  in  the  gas-surface  interaction  that  can  only  be 
understood  within  a  quantum  framework,  namely  elastic  scattering  and  diffraction. 

There  are  two  ways  to  determine  quantum  mechanical  reaction  probabilities:  by  solving  the  time- 
dependent  or  the  time- in  dependent  Schrodinger  equation.  Both  approaches  are  equivalent  [9]  and 
should  give  the  same  results.  The  answer  to  the  question  which  method  is  more  appropriate  depends 
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Figure  2:  Illustration  of  the  different  processes  that  can  occur  when  atoms  or  molecules  are  imping¬ 
ing  on  a  surface.  The  substrate  with  lattice  constant  a  is  represented  in  a  simple  ball  and  spring 
picture. 


on  the  particular  problem.  Time-independent  implementations  are  usually  more  restrictive  as  far 
as  the  form  of  the  potential  is  concerned,  but  often  the  choice  of  the  method  is  a  matter  of  training 
and  personal  taste. 


In  the  time-dependent  or  wave-packet  formulation,  the  solution  of  the  time-dependent  Schrodinger 
equation 

ih^*(R,t)=H*(R,t)  (11) 

can  formally  be  written  as 

\k(R,  t)  =  e-iHt/h  (R,  t  =  0)  ,  (12) 

if  the  potential  is  time-independent.  The  most  common  methods  to  represent  the  time-evolution 
operator  exp {—iHt/K)  in  the  gas-surface  dynamics  community  are  the  split-operator  [10,  11]  and 
the  Chebychev  [12]  methods.  In  the  split-operator  method,  the  time-evolution  operator  for  small 
time  steps  At  is  written  as 


e-iHAt/h  _  e-iKAt/2h  e~iVAt/h 


e-iKAt/2h  +  0(At3)  ^ 


(13) 


where  K  is  the  kinetic  energy  operator  and  V  the  potential  term.  In  an  alternative  approach,  the 
Chebyshev  method,  the  time-evolution  operator  is  expanded  as 

Jmax 

e-,HA,/H  =  a, (At)  Tj(H)  , 

3= 1 

where  the  Tj  are  Chebyshev  polynomials  and  H  is  the  Hamiltonian  rescaled  to  have  eigenvalues 
in  the  range  (—1,1).  Both  propagation  schemes  use  the  fact  that  the  kinetic  energy  operator 
is  diagonal  in  k-space  and  the  potential  is  diagonal  in  real-space.  The  wave  function  and  the 
potential  are  represented  on  a  numerical  grid,  and  the  switching  between  the  k-space  and  real- 
space  representations  is  efficiently  done  by  Fast  Fourier  Transformations  (FFT)  [5]. 


(14) 
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In  the  time-independent  formulation,  on  the  other  hand,  the  wave  function  is  usually  expanded 
in  a  suitable  basis  set.  This  often  requires  to  introduce  the  concept  of  reaction  path  coordinates. 
Starting  from  the  time-independent  Schrodinger  equation 

(H-E)'&  =  0,  (15) 

one  specific  reaction  path  coordinate  s  is  chosen.  Then  the  kinetic  energy  operator  in  this  coordinate 
is  separated  yielding 

(^d2s+HE)*  =  0  .  (16) 

Here  H  is  the  original  Hamiltonian  except  for  the  kinetic  energy  operator  in  the  reaction  path 
coordinate.  Usually  the  use  of  curve-linear  reaction  path  coordinates  results  in  a  more  complicated 
expression  for  the  kinetic  energy  operator  involving  cross  terms,  but  for  the  sake  of  clarity  this  has 
been  neglected  in  Eq.  (16).  In  the  coordinates  perpendicular  to  the  reaction  path  coordinate,  the 
wave  function  is  expanded  in  some  suitable  set  of  basis  functions, 

tf  =  tf(s,...)  =  Ipn{s)  \n)  .  (17) 

n 

Here  n  is  a  multi-index,  and  the  expansion  coefficients  il>n(s)  are  assumed  to  be  a  function  of  the 
reaction  path  coordinate.  This  expansion  of  is  inserted  in  the  Schrodinger  equation  (16),  and  this 
equation  is  multiplied  by  {m\  which  corresponds  to  performing  a  multi-dimensional  integral.  Since 
the  basis  functions  |  n)  are  assumed  to  be  independent  of  s,  this  yields  the  so-called  coupled-channel 
equations 

Y]  |  -  E)  5m,n  +  (m\H\n)  \  i>n(s)  =  0  .  (18) 

Instead  of  a  high-dimensional  partial  differential  equation  -  the  original  time- independent  Schrodin- 
ger  equation  (15)  -  there  is  now  a  set  of  coupled  ordinary  differential  equation.  Still  a  straight¬ 
forward  numerical  integration  of  the  coupled-channel  equations  leads  to  instabilities,  except  for 
in  simple  cases,  due  to  exponentially  increasing  so-called  closed  channels.  These  problems  can  be 
avoided,  for  example  by  making  the  potential  stepwise  constant  so  that  the  wave  function  can  be 
analytically  propagated  [13,  14] 


3  Scattering  at  surfaces 


In  Fig.  2,  possible  collision  processes  in  the  scattering  of  atoms  and  molecules  at  surfaces  are 
summarized.  A  monoenergetic  beam  of  atoms  or  molecules  characterized  by  the  wave  vector  K*  = 
P i/h  where  P,  is  the  initial  momentum  of  the  particles  is  impinging  on  a  periodic  surface  with 
lattice  constant  a.  Classically  there  will  always  be  a  certain  energy  transfer  from  the  molecules 
to  the  substrate  when  the  incoming  particles  hit  the  surface.  Quantum  mechanically,  however, 
there  will  be  a  non-vanishing  probability  for  elastic  scattering,  i.e.  with  no  energy  transfer  to  the 
substrate.  This  probability  is  given  by  the  so-called  Debye- Waller  factor. 

Furthermore,  if  the  de  Broglie  wavelength  A  =  l/|Kj|  of  the  incident  beam  is  of  the  order  of 
the  lattice  spacing  a,  quantum  effects  in  the  momentum  transfer  parallel  to  the  surface  become 
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important.  The  periodicity  of  the  substrate  leads  to  the  conservation  of  the  quasi-momentum 
parallel  to  the  surface  which  means  that  the  component  of  the  wave  vector  parallel  to  the  surface 
can  only  be  changed  by  a  reciprocal  lattice  vector  of  the  periodic  surface: 

Kj  =  K|  +  Gmn  ,  (19) 

where  Gmn  is  a  reciprocal  lattice  vector  of  the  periodic  surface  and  k|  and  kJ|  are  the  initial 
and  final  wave  vectors  parallel  to  the  surface.  As  a  consequence,  diffraction  results,  i.e. ,  there  is 
only  a  discrete  number  of  allowed  scattering  angles.  The  intensity  of  the  elastic  diffraction  peak 
mn  according  to  Eq.  (19)  is  denoted  by  Imn.  The  scattering  peak  loo  with  kJ|  =  k|  is  called  the 
specular  peak.  It  is  important  to  note  that  the  allowed  scattering  angles  are  a  consequence  of  the 
surface  geometry  which  means  that  from  the  diffraction  pattern  the  periodicity  and  lattice  constant 
of  the  surface  can  be  derived,  whereas  the  intensity  of  the  different  diffraction  peaks  depends  on 
the  particular  molecule-surface  system.  Therefore  the  coherent  scattering  of  atoms  or  molecules 
from  surfaces  can  be  used  as  a  tool  for  probing  surface  structures  which  has  first  been  realized  in 
1930  [15].  In  particular  helium  atom  scattering  (HAS)  experiments  have  been  carried  out  to  study 
the  surface  crystallography  (see,  e.g.,  [16,  17]  and  references  therein). 

The  probability  for  elastic  scattering,  i.e.  the  Debye- Waller  factor,  vanishes  rapidly  with  increasing 
mass  of  the  impinging  molecules.  For  heavier  molecules,  predominantly  inelastic  scattering  occurs. 
The  main  source  for  the  energy  transfer  between  the  impinging  molecules  and  the  substrate  is  the 
excitation  and  deexcitation  of  substrate  phonons  but  also  electron-hole  pairs  of  the  substrate  may 
be  involved.  Since  phonons  also  carry  momentum,  the  conservation  of  quasi-momentum  parallel 
to  the  surface  is  modified  to 


K1}  =  K|  +  Gmn  +  Y  ±Q  > 

exch.phon. 


(20) 


where  Q  is  a  two-dimensional  phonon-momentum  vector  parallel  to  the  surface.  The  plus-signs  in 
the  sum  correspond  to  the  excitation  or  emission  of  a  phonon  while  the  minus-signs  represent  the 
deexcitation  or  absorption  of  a  phonon.  The  energy  balance  in  phonon-inelastic  scattering  can  be 
expressed  as 


2  M 


ft2K? 
2  M 


Y  ±hujQJ  ■ 


(21) 


exch.phon. 


where  huo qj  corresponds  to  the  energy  of  the  phonon  with  momentum  Q  and  mode  index  j . 
From  the  change  of  the  momentum  and  the  energy  in  single-phonon  scattering,  the  surface  phonon 
spectrum  can  be  derived  which  has  been  extensively  done  using  helium  atom  scattering  as  a 
probe  [16,  18].  The  excitation  of  phonons  usually  leads  to  a  reduced  normal  component  of  the 
kinetic  energy  of  the  back-scattered  atoms  or  molecules.  Thus  the  reflected  beam  is  shifted  in 
general  to  larger  angles  with  respect  to  the  surface  normal  compared  to  the  angle  of  incidence. 
The  resulting  supraspecular  scattering  is  indicated  in  Fig.  2  as  the  inelastic  reflection  event. 


In  addition  to  diffraction,  often  resonances  in  the  intensity  of  the  specular  peak  as  a  function  of 
the  angle  of  incidence  are  observed  in  the  case  of  the  scattering  of  weakly  interacting  particles  at 
smooth  surfaces,  [19].  These  so-called  selective  adsorption  resonances  are  also  indicated  in  Fig.  2. 
They  occur  when  the  scattered  particle  can  make  a  transition  into  one  of  the  discrete  bound  state 
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of  the  adsorption  potential  [20].  This  is  only  possible  if  temporarily  the  motion  of  the  particle 
is  entirely  parallel  to  the  surface.  The  interference  of  different  possible  paths  along  the  surface 
causes  the  resonance  effects.  Energy  and  momentum  conservation  yields  the  selective  adsorption 
condition 

h2Kj  =  n2(  Kj  +  Gmn)2 
2.1/  2.1/ 

where  Ei  is  a  bound  level  of  the  adsorption  potential.  From  the  scattering  resonances,  bound 
state  energies  can  be  obtained  using  Eq.  (22)  without  any  detailed  knowledge  about  the  scattering 
process. 


Since  molecules,  in  contrast  to  atoms,  have  also  internal  degrees  of  freedoms,  namely  rotations 
and  vibrations,  additional  peaks  may  appear  in  the  diffraction  pattern  of  molecules.  They  are  due 
to  the  fact  that  these  rotations  and  vibrations  can  be  excited  and  de-excited  during  the  collision 
process.  The  total  energy  balance  in  the  molecular  scattering  is  therefore  more  complex: 


2  M 


n2K2 

2  M 


+  AErot  +  AEvib 


+  ^2  ±^Qj  ■ 

exch.phon. 


(23) 


Since  the  time-scale  of  the  molecular  vibrations  is  usually  much  shorter  than  the  scattering  time 
or  the  rotational  period  and  therefore  the  associated  excitation  energies  are  much  larger,  the 
excitation  of  molecular  vibrations  in  molecule-surface  scattering  is  usually  negligible,  in  contrast 
to  the  phonon  excitation.  Molecular  rotations,  on  the  other  hand,  can  be  excited  rather  efficiently 
in  the  scattering  at  highly  corrugated  and  anisotropic  surfaces.  This  leads  to  additional  peaks  in 
the  diffraction  spectrum,  the  rotationally  inelastic  diffraction  peaks. 


Experimentally,  rotationally  inelastic  diffraction  of  hydrogen  molecules  has  been  first  observed  in 
the  scattering  at  inert  ionic  solids  such  as  MgO  [21]  or  NaF  [22],  At  metal  surfaces  with  a  high 
barrier  for  dissociative  adsorption,  the  molecules  are  scattered  at  the  tails  of  the  metal  electron 
density  which  are  usually  rather  smooth.  Hence  relatively  weak  diffraction  and  hardly  any  rota¬ 
tionally  inelastic  diffraction  has  been  observed  for,  e.g.,  the  scattering  of  H2  from  Cu(001)  [23,  24]. 
This  is  different  for  the  case  of  HD  scattering,  where  the  displacement  of  the  center  of  mass  from 
the  center  of  the  charge  distribution  leads  to  a  strong  rotational  anisotropy  [25]. 

At  reactive  surfaces  where  non-activated  adsorption  is  possible,  the  scattering  occurs  rather  close 
to  the  surface  where  the  potential  energy  surface  is  already  strongly  corrugated  and  anisotropic. 
For  such  systems,  intensive  rotationally  inelastic  peaks  in  the  diffraction  pattern  have  been  pre¬ 
dicted  theoretically  in  six- dimensional  quantum  dynamical  calculations  [26]  for  the  scattering  of 
H2/Pd(100)  where  non-activated  together  with  activated  pathways  to  dissociative  adsorption  ex¬ 
ist.  One  typical  calculated  angular  distribution  of  H2  molecules  scattered  at  Pd(100)  is  shown 
in  Fig.  3  [26].  The  total  initial  kinetic  energy  is  Ei  =  76  meV.  The  incident  parallel  momentum 
equals  2 hG  along  the  (Oil)  direction  which  corresponds  to  an  incident  angle  of  6j  =  32°.  The 
molecules  are  initially  in  the  rotational  ground  state  jt  =  0.  Figure  3a  shows  the  so-called  in-plane 
scattering  distribution,  i.e.  the  diffraction  peaks  in  the  plane  spanned  by  the  wave  vector  of  the 
incident  beam  and  the  surface  normal.  The  label  (m,  n)  denotes  the  parallel  momentum  transfer 
AG ||  =  ( mG,nG ).  The  specular  peak  is  the  most  pronounced  one,  but  the  first  order  diffraction 
peak  (10)  is  only  a  factor  of  four  smaller.  Note  that  in  a  typical  helium  atom  scattering  experiment 
the  off-specular  peaks  are  about  two  orders  smaller  than  the  specular  peak  [16].  This  is  due  to 
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Figure  3:  Diffraction  spectrum  of  H2  scattered  at  Pd(100)  for  a  kinetic  energy  of  76  meV  at 
an  incidence  angle  of  32°  along  the  [10]  direction  of  the  square  surface  lattice,  obtained  by  six¬ 
dimensional  quantum  coupled  channel  calculations,  a)  In-plane  diffraction  spectrum  where  all  peaks 
have  been  labeled  according  to  the  transition,  b)  In-plane  and  out-of-plane  diffraction  peaks.  The 
open  and  filled  circles  correspond  to  the  rotationally  elastic  and  rotationally  inelastic  scattering, 
respectively,  with  he  radius  of  the  circles  being  proportional  to  the  logarithm  of  the  scattering 
intensity  (after  [26]). 


the  fact  that  the  chemically  inert  helium  atoms  are  scattered  at  the  smooth  tails  of  the  surface 
electron  distribution. 

The  peaks  labeled  with  A j  =  2  in  Fig.  3a  correspond  to  rotationally  inelastic  diffraction  involving 
the  rotational  excitation  j  =  0  — >  2  summed  up  over  all  final  azimuthal  quantum  numbers  rrij. 
Except  for  the  specular  peak,  the  intensities  of  the  rotationally  inelastic  diffraction  peaks  are 
even  larger  than  the  corresponding  rotationally  elastic  diffraction  peaks  with  the  same  momentum 
transfer  ( m ,  n).  Note  that  because  of  the  particular  conditions  with  the  incident  parallel  momentum 
corresponding  to  the  reciprocal  lattice  vector  Gy  =  (2G,  0),  the  rotationally  elastic  and  inelastic 
(20)  diffraction  peaks  fall  upon  each  other. 

The  full  diffraction  pattern  including  also  the  so-called  out-of-plane  scattering  peaks  is  shown 
in  Fig.  3b.  The  open  circles  represent  the  rotationally  elastic,  the  filled  circles  the  rotationally 
inelastic  diffraction  peaks.  The  radii  of  the  circles  are  proportional  to  the  logarithm  of  the  scattering 
intensity.  It  is  obvious  that  there  is  a  significant  fraction  of  out-of-plane  scattering  with  the  sum 
of  all  out-of-plane  scattering  intensities  approximately  equal  to  the  sum  of  all  in-plane  scattering 
intensities.  Interestingly,  some  diffraction  peaks  with  a  large  parallel  momentum  transfer  still  show 
substantial  intensities.  This  phenomenon  is  well  known  from  helium  atom  scattering  and  has  been 
discussed  within  the  concept  of  so-called  rainbow  scattering  [27] . 

Experimentally,  it  is  not  so  easy  to  measure  diffraction  patterns  at  reactive  surfaces  since  during  the 
experiment,  a  significant  fraction  of  the  impinging  molecules  remains  on  the  surface  which  destroys 
the  ideal  periodicity  of  the  surface  and  therefore  also  suppress  the  occurrence  of  diffraction  patterns. 
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Thus  the  surface  has  to  be  kept  clean  which  is  achieved  by  relatively  high  surface  temperatures  so 
that  adsorbates  quickly  desorb  again.  High  surface  temperatures,  on  the  other  hand,  also  smear 
out  the  diffraction  pattern.  Still  rotationally  inelastic  peaks  in  addition  to  rotationally  elastic  peaks 
have  been  clearly  identified  in  the  diffraction  pattern  of  D2/Ni(110)  [28]  and  D2/Rh(110)  [29]. 

At  reactive  surfaces,  the  particles  can  off  course  also  adsorb.  As  it  is  indicated  in  Fig.  2,  molecules 
can  adsorb  both  molecularly  which  means  intact  or  dissociatively.  In  the  case  of  the  atomic  or 
molecular  adsorption,  the  particles  can  only  remain  trapped  at  the  surface  if  their  initial  kinetic  is 
transfered  to  the  surface  and  dissipated.  For  light  projectiles,  the  quantum  nature  of  the  substrate 
phonons  becomes  important  in  the  energy  transfer  process.  These  topics  will  be  discussed  in  the 
next  section. 

4  Atomic  and  molecular  adsorption 


The  sticking  or  adsorption  probability  is  defined  as  the  fraction  of  atoms  or  molecules  impinging  on 
a  surface  that  are  not  scattered  back,  i.e.  that  remain  on  the  surface.  It  is  important  to  note  that  in 
the  case  of  atomic  or  molecular  adsorption  when  the  molecule  stays  intact,  the  particles  can  only 
remain  on  the  surface  if  they  transfer  their  energy  to  the  substrate.  This  is  similar  to  gas-phase 
reactions  where  a  bond  between  two  reactants  can  only  be  formed  in  a  three-body  collision  where 
a  third  reaction  partner  has  to  carry  away  the  energy  gained  by  the  reaction  unless  there  are  other 
dissipation  channels  such  as  radiation. 

At  a  surface,  there  are  two  main  channels  for  energy  dissipation  namely  phonon  and  electron-hole 
pair  excitations.  Here  we  focus  on  the  energy  transfer  to  phonons  since  they  usually  represent 
the  main  channel  for  dissipation  [30].  For  the  explicit  evaluation  of  sticking  probabilities,  Pg(e)  is 
defined  as  the  probability  that  an  incoming  particle  with  kinetic  energy  E  transfers  the  energy  e 
to  the  surface.  Only  if  the  particle  transfers  more  than  its  initial  energy  to  substrate  excitations, 
it  can  remain  at  the  surface.  Hence  the  sticking  probability  can  be  expressed  as 

/»oo 

S(E)  =  /  PE(e)  de.  (24) 

J  E 


In  order  to  discuss  the  essentials  of  atomic  and  molecular  adsorption  due  to  the  energy  transfer 
to  phonons,  we  will  use  the  rather  simple  hard-cube  model  (HCM)  [31,  32],  In  this  model  that  is 
schematically  illustrated  in  Fig.  4a,  the  impact  of  the  atom  on  the  surface  is  treated  as  a  binary 
elastic  collision  between  a  gas  phase  atom  (mass  rn)  and  a  substrate  atom  (mass  Mc)  which  is 
moving  freely  with  a  velocity  distribution  Pc(vc).  Because  of  the  adsorption  well  of  depth  E^d,  the 
particle  is  accelerated  and  impinges  on  the  hard  cube  with  a  velocity 


^well  =  ~\  K  + 


2  Fad 


m 


(25) 


Due  to  the  simplicity  of  the  hard  cube  model,  it  can  be  solved  analytically  [32].  Assuming  a 
weighted  Maxwellian  velocity  distribution  for  vc,  the  trapping  probability  becomes 


Strap(vg)  =  \  +  \  erf(auiim)  + 


expj-o^L} 

2 \Je  o,/cwe]] 


(26) 
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Figure  4:  Hard  cube  model,  a)  Schematic  illustration  of  the  model.  An  atom  or  molecule  with 
mass  m  is  impinging  in  an  attractive  potential  with  well  depth  Fad  on  a  surface  modeled  by  a 
cube  of  effective  mass  Mc.  The  surface  cube  is  moving  with  a  velocity  vc  given  by  a  Maxwellian 
distribution,  b)  Trapping  probability  as  a  function  of  the  kinetic  energy  evaluated  according  to  the 
hard  cube  model  Eq.  (26)  for  different  adsorption  energies  Fa a,  mass  ratios  fi  =  m/Mc  and  surface 
temperatures  Ts. 


where  a  =  Mc/2ksTs ,  viim  is  given  by 

_  /i  + 1  /2Ead  /x  -  1 

Vlim  ~  2  V  m  2  Vwe11  ’ 

and  /i  is  the  mass  ratio  /x  =  m/M. 


(27) 


Typical  sticking  curves  obtained  with  the  hard  cube  model  are  shown  in  Fig.  4b.  All  curves  show 
the  same  behavior,  namely  a  monotonic  decrease  of  the  sticking  probability  as  a  function  of  the 
kinetic  energy.  This  dependence  is  a  consequence  of  the  fact  that  the  energy  transfer  to  the  surface 
becomes  less  efficient  at  higher  kinetic  energies.  Although  more  energy  is  transfered  to  the  surface 
at  higher  kinetic  energies,  the  fraction  of  particles  that  loose  more  energy  than  their  initial  kinetic 
energy  becomes  smaller. 


The  examples  shown  in  Fig.  4b  have  been  chosen  in  order  to  illustrate  further  general  trends  in 
atomic  adsorption.  If  the  mass  ratio  between  the  impinging  molecules  and  the  substrate  atoms 
increases,  i.e.  if  heavier  atoms  are  hitting  the  surface,  more  energy  is  transfered  to  the  surface 
so  that  the  trapping  probability  increases.  If  the  adsorption  well  becomes  deeper,  the  impinging 
molecules  become  faster  and  loose  more  energy  upon  impact  which  also  increases  the  trapping 
probability.  For  an  adsorption  well  of  0.6  eV,  the  effect  of  changing  the  surface  temperature  is 
shown.  A  higher  surface  temperature  Ts  leads  to  a  broader  velocity  distribution  which  results  in 
an  averaging  over  a  wider  range  of  kinetic  energies.  This  causes  a  decrease  for  negative  curvature 
of  the  sticking  curve,  i.e.  at  high  sticking  probabilities,  and  an  increase  for  positive  curvature,  i.e. 
at  low  sticking  probabilities. 


Adsorption  well  depths  have  often  been  estimated  by  fitting  measured  trapping  probabilities  to 
curves  derived  from  the  hard  cube  model,  also  for  molecular  adsorption  (see,  e.g.,  [32,  33]).  In  the 
simple  hard  cube  model  the  surface  is  flat  and  structureless  which  means  that  in  any  scattering 
and  adsorption  process  the  incident  parallel  momentum  is  conserved.  For  the  sticking  probability 
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Figure  5:  Sticking  probability  of  rare  gas  atoms  on  Ru(0001)  at  a  surface  temperature  of  Ts  =  6.5  K. 
Stars  (*):  experiment;  lines:  theoretical  results  obtained  with  the  forced  oscillator  model  (after  [34], 
not  all  measured  data  points  are  included) 


this  leads  to  the  normal  energy  scaling ,  i.e. ,  the  sticking  probability  is  a  function  of  the  normal 
component  E{  cos 2  6i  of  the  incident  energy  alone,  where  6i  is  the  angle  of  incidence.  Real  surfaces, 
however,  are  not  structureless  as  far  as  the  interaction  of  atoms  and  molecules  is  concerned  since 
adsorption  corresponds  to  the  making  of  a  chemical  bond  which  strongly  depends  on  the  local 
environment.  This  leads  to  corrugation  in  the  potential  energy  surface,  i.e.,  the  potential  depends 
on  the  lateral  position  of  the  interacting  particle  on  the  surface.  Trapping  probabilities  often  scale 
as  Ei  cosn  6i  with  n  <  2.  An  exponent  of  n  =  0  corresponds  to  total  energy  scaling which  is  usually 
associated  with  a  highly  corrugated  potential  energy  surface. 

A  further  prediction  of  the  hard  cube  model  is  that  the  sticking  probability  in  the  limit  E  — >  0  and 
Ts  — >  0  always  becomes  unity  no  matter  how  small  the  adsorption  well,  no  matter  how  small  the 
mass  ratio  between  the  impinging  atom  and  the  substrate  oscillator,  if  there  is  no  barrier  before 
the  adsorption  well.  This  is  due  to  the  energy  transfer  to  the  surface  which  in  the  limit  of  zero 
energy  does  not  allow  the  particle  to  return  into  the  gas-phase. 

Experiments  indeed  show  that  for  many  adsorption  systems  all  impinging  particles  become  trapped 
at  low  kinetic  energies,  for  example  in  the  case  of  Xe  impinging  on  Ru(0001)  [35,  34]  which  is  shown 
in  Fig.  5.  However,  for  the  lighter  rare  gases  Kr,  Ar  and  in  particular  Ne  the  trapping  probability 
extrapolated  to  E  — *  0  is  clearly  below  one.  This  behavior  can  not  be  reproduced  using  classical 
mechanics  [35] ,  instead  it  can  only  be  understood  if  the  quantum  nature  of  the  substrate  phonon 
system  is  taken  into  account  leading  to  a  non-zero  probability  for  elastic  scattering  at  the  surface. 
A  classical  treatment  of  the  solid  is  only  appropriate  if  the  energy  transfer  to  the  surface  is  large 
compared  to  the  Debye  energy  of  the  solid  [36] . 

The  essential  physics  can  be  captured  by  just  considering  an  atomic  projectile  interacting  via 
linear  coupling  with  a  single  quantum  surface  oscillator  in  the  so-called  trajectory  approximation. 
In  this  model,  the  classical  motion  of  the  incoming  particle  is  assumed  not  to  be  perturbed  by  the 
coupling  to  the  surface  oscillator.  The  classical  trajectory  then  introduces  a  time-dependent  force 
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in  the  Hamiltonian  of  the  oscillator 

Hosc  =  fru)  (a+o  +  — )  +  A (t)  (a+  +  a)  ,  (28) 

In  this  forced  oscillator  model  [37,  38,  39],  the  mean  number  of  excited  phonon  n  is  given  basically 
by  the  square  of  the  Fourier  transform  of  the  coupling  A  (t), 


n 


\  J  e-  m# 


(29) 


The  number  n  also  enters  the  expression  of  probability  P]t  for  a  transition  from  an  initial  oscillator 
state  i  to  the  final  state  j  [37]: 

Pji  =  j,  e~n  ni-*  [L]~\n)}2  ,  j>i.  (30) 

Here  L3-  1  is  an  associated  Laguerre  polynomial  [40].  For  an  excitation  PjO  from  the  ground  state 
one  just  obtains  a  Poisson  distribution.  The  probabilities  (30)  yield  the  energy  distribution  of 
excited  phonons  and  hence  also  the  energy  transfer  from  the  impinging  particle  to  the  surface 
which  that  they  correspond  to  the  probability  Pe{ e)  entering  the  expression  (24)  for  the  sticking 
probability. 

This  forced  oscillator  model  has  been  applied  to  evaluate  the  sticking  probabilities  of  rare  gas 
atoms  on  a  Ru(001)  surface  at  a  temperature  of  Ts  =  6.5  K  [34],  Using  the  potential  well  depth, 
the  potential  range,  the  mass  of  the  surface  oscillator  and  the  surface  Debye  temperature  as 
parameters,  this  simple  model  was  able  reproduce  the  measured  data  quite  well  (see  Fig.  5).  In 
particular  for  the  a  light  rare  gas  such  as  Ne  impinging  on  a  metal  surface,  the  quantum  effects  in 
the  surface  recoil  are  quite  substantial. 

At  even  lower  kinetic  energies  than  reached  in  the  experiments  [34]  shown  in  Fig.  5,  the  quantum 
nature  of  the  impinging  particles  cannot  be  neglected  any  longer.  Hence  the  trajectory  approxi¬ 
mation  cannot  be  applied  any  more.  In  fact,  in  the  limit  E  — »  0  the  de  Broglie  wavelength  of  the 
incoming  particle  tends  to  infinity.  In  the  case  of  a  short-range  attractive  potential  this  means  that 
the  amplitude  of  the  particle’s  wave  function  vanishes  in  the  attractive  region  [36,  41].  Thus  there  is 
no  coupling  and  consequently  no  energy  transfer  between  the  particle  and  the  substrate  vibrations. 
Therefore  the  quantum  mechanical  sticking  probability  also  vanishes  for  E  — »  0.  However,  in  order 
to  see  this  effect  extremely  small  kinetic  energies  corresponding  to  a  temperature  below  0.1  K  are 
required  [36].  Nevertheless,  this  quantum  phenomenon  in  the  sticking  at  surfaces  has  been  verified 
experimentally  for  the  adsorption  of  atomic  hydrogen  on  thick  liquid  4He  films  [42]. 


5  Dissociative  adsorption 


In  the  case  of  dissociative  adsorption  there  is  another  channel  for  energy  transfer,  namely  the 
conversion  of  the  kinetic  and  internal  energy  of  the  molecule  into  translational  energy  of  the  atomic 
fragments  on  the  surface  relative  to  each  other.  This  bond- breaking  process  which  is  illustrated  in 
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Figure  6:  Contour  plots  of  the  potential  energy  surface  along  two-dimensional  cuts  through  the 
six-dimensional  coordinate  space  of  H2  in  front  of  (100)  metal  surfaces  determined  by  DFT-GGA 
calculations  for  H2  molecules  above  the  bridge  site.  The  contour  spacing  is  0.1  eV  per  H2  molecule, 
(a)  H2/Pd(100)  (after  [43]),  (b)  H2/Cu(100)  [44],  The  inset  is  a)  illustrates  the  process  of  disso¬ 
ciative  adsorption  while  in  b)  two  typical  trajectories  showing  the  effect  of  initial  vibrations  are 
included. 


the  inset  of  Fig.  6  represents  the  fundamental  difference  to  atomic  or  molecular  adsorption.  As  far 
as  the  complete  adsorption  process  is  concerned,  it  is  true  that  eventually  the  atomic  fragments 
will  also  dissipate  their  kinetic  energy  and  come  to  rest  at  the  surface.  However,  in  the  case  of  light 
molecules,  in  particular  molecular  hydrogen,  dissociating  on  metal  surfaces  the  energy  transfer 
to  the  substrate  is  very  small  because  of  the  large  mass  mismatch.  Hence  the  probability  for 
dissociative  adsorption  is  almost  entirely  given  by  the  initial  dissociation  probability  upon  the 
impact  on  the  surface  since  the  fragments  do  in  general  not  directly  recombine  and  desorb  again. 
Therefore  the  dissociative  adsorption  process  can  be  described  within  low-dinrensional  potential 
energy  surfaces  neglecting  the  surface  degrees  of  freedom.  This  also  requires  that  there  is  no 
substantial  surface  rearrangement  upon  adsorption,  but  this  is  usually  the  case  in  the  dissociative 
adsorption  on  close-packed  metal  surfaces. 

The  dynamics  of  the  interaction  of  hydrogen  with  metal  surfaces  has  been  well-studied,  both  ex¬ 
perimentally  [45]  and  theoretically  [2,  46,  47,  48].  In  Fig.  6,  two  cuts  through  the  six-dimensional 
potential  energy  surface  derived  from  DFT  calculations  for  two  benchmark  systems  are  plotted, 
for  H2/Pd(100)  where  non-activated  dissociation  is  possible  [43])  and  for  the  activated  system 
H2/Cu(100)  [44],  The  cuts  including  the  respective  minimum  energy  paths  to  dissociative  adsorp¬ 
tion  correspond  to  so-called  elbow  plots  which  show  the  potential  energy  as  a  function  of  the  H2 
center  of  mass  distance  from  the  surface  and  the  H-H  distance. 

In  order  to  quantitatively  evaluate  the  adsorption  probability  of  hydrogen,  the  full  six-dimensional 
PES  has  to  be  known.  Furthermore,  the  dynamics  of  hydrogen  require  a  quantum  treatment  be¬ 
cause  of  its  light  mass.  Because  of  the  high  computational  effort  associated  with  quantum  methods, 
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Figure  7:  Dissociative  adsorption  of  H2/Pd(100).  a)  Sticking  probability  of  H2/Pd(100)  as  a  func¬ 
tion  of  the  initial  kinetic  energy.  Circles:  experiment  [51],  dashed  and  solid  line:  theory  according 
to  H2  initially  in  the  ground  state  and  with  a  thermal  distribution  appropriate  for  a  molecular 
beam  [49].  The  inset  shows  the  theoretical  results  using  an  improved  ab  initio  potential  energy 
surface  [52],  b)  Illustration  of  the  steering  effect. 


for  a  long  time  the  theoretical  treatment  was  limited  to  studies  within  a  reduced  dimensionality. 
Only  recently  the  first  quantum  studies  were  performed  in  which  the  full  dimensionality  of  the 
hydrogen  molecule  was  taken  into  account  [49,  50]. 

The  corresponding  dissociative  adsorption  probability  of  H2/Pd(100)  is  shown  in  Fig.  7a.  Experi¬ 
ment  [51]  and  theory  [49,  52]  agree  well,  as  far  as  the  qualitative  trend  of  the  adsorption  probability 
as  a  function  of  the  kinetic  energy  is  concerned.  First  there  is  an  initial  decrease,  and  after  a  min¬ 
imum  the  sticking  probability  rises  again.  The  initial  decrease  of  the  sticking  probability  is  typical 
for  H2  adsorption  at  transition  metal  surfaces  [51,  53,  45,  54,  55].  Originally,  such  an  initial  de¬ 
crease  was  associated  with  the  existence  of  a  molecular  adsorption  precursor  state  through  which 
the  dissociation  was  assumed  to  proceed.  However,  nowadays  it  is  well  established  that  this  be¬ 
havior  in  the  dissociation  of  H2  on  reactive  metal  surfaces  is  caused  by  dynamical  effects,  namely 
steering  [56,  49,  57]  and  dynamical  trapping  [58,  59,  60]. 

The  steering  effect  is  illustrated  in  Fig.  7b  where  typical  trajectories  are  plotted  as  a  function  of 
a  lateral  coordinate  and  a  reaction  path  coordinate  which  connects  the  molecule  in  the  gas  phase 
with  the  dissociated  molecule  on  the  surface.  It  is  important  to  note  that  in  these  systems  the 
PES  shows  purely  attractive  paths  towards  dissociative  adsorption,  but  the  majority  of  reaction 
paths  for  different  molecular  orientations  and  impact  points  exhibits  energetic  barriers  hindering 
the  dissociation.  Still,  at  low  kinetic  energies  most  impinging  molecules  are  steered  towards  the 
attractive  dissociation  channel  leading  to  a  high  adsorption  probability.  In  addition,  the  corrugation 
and  anisotropy  leads  to  a  conversion  of  the  kinetic  energy  of  the  impinging  molecule  into  internal 
molecular  degrees  of  freedom  such  as  rotation  and  vibration  and  into  lateral  motion  along  the 
surface.  This  energy  is  then  not  available  for  a  direct  escape  from  the  adsorption  well  so  that 
the  molecules  become  dynamically  trapped.  Both  mechanisms  that  are  only  operative  a  low  kinetic 
energies  result  large  adsorption  probabilities,  and  their  suppression  for  higher  kinetic  energies  leads 
to  the  decrease  of  the  sticking  probability.  At  even  higher  energies,  the  molecules  can  directly  cross 
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Figure  8:  Dissociative  adsorption  probability  of  H2  on  Cu(100)  as  a  function  of  the  incident  kinetic 
energy  determined  by  six-dinrensional  quantum  wave-packet  calculations  for  molecules  initially  in 
the  vibrational  ground  state  and  first  excited  state,  respectively  [61].  The  calculations  are  compared 
to  experimental  results  derived  from  an  analysis  of  adsorption  and  desorption  experiments  [62]. 


the  dissociation  barriers  causing  again  an  increase  in  the  adsorption  probability. 

As  far  as  the  activated  adsorption  is  concerned,  the  interaction  of  hydrogen  with  copper  surfaces 
has  served  as  a  model  system  [46,  47,  63,  64,  50,  65,  66].  It  was  also  the  first  system  for  which 
high-dinrensional  potential  energy  surface  were  mapped  out  by  DFT  methods  [67,  68,  69].  Figure  8 
shows  the  calculated  sticking  probabilities  for  molecules  either  in  the  vibrational  ground  state  or 
in  the  first  excited  state,  respectively,  based  on  an  six-dinrensional  wave-packet  calculations  on  an 
ab  initio  PES  [50,  61].  The  theoretical  results  are  rather  close  to  the  experimental  data  which  were 
derived  from  an  analysis  of  both  adsorption  and  desorption  experiments  [62]  by  using  the  following 
analytical  form  of  the  vibrationally  resolved  sticking  probability  as  a  function  of  the  kinetic  energy: 

S„(£)  =  4{1+tanh(t|^))}  (31) 

The  agreement  between  theory  and  experiment  in  Fig.  8  is  very  satisfactory  except  for  the  high- 
energy  behavior.  However,  the  experimental  data  of  H2  in  Fig.  8  for  kinetic  energies  above  0.5  eV 
were  derived  from  thermal  desorption  experiments.  As  far  the  saturation  value  of  the  experimental 
sticking  probability  is  concerned,  there  is  a  large  uncertainty  since  the  high  energy  contributions 
in  desorption  are  exponentially  suppressed  through  the  Boltzmann  factor. 

The  sticking  probability  of  H2/Cu(100)  shows  a  typical  activated  behavior  with  an  onset  at  ap¬ 
proximately  0.5  eV  for  H2  for  molecules  initially  in  the  vibrational  ground  state.  This  onset  is  given 
by  the  minimum  energy  barrier  including  zero-point  effects  which  arise  from  the  quantization  of 
the  molecular  levels  due  to  the  localization  of  the  wave  function  at  the  minimum  barrier  position. 
The  slope  of  the  sticking  probability  is  directly  related  to  the  distribution  of  the  barrier  heights 
for  dissociative  adsorption  in  the  multidimensional  potential  energy  surface  [64] .  Thus  sticking  can 
be  understood  in  terms  of  the  region  of  the  surface  that  classically  is  available  to  dissociation. 
This  so-called  hole  model  [70]  is  valid  at  high  kinetic  energies  when  the  incoming  particles  are  not 
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significantly  redirected  by  the  shape  of  the  potential  energy  surface. 


As  Fig.  8  demonstrates,  for  initially  vibrationally  excited  molecules  the  sticking  probability  is 
significantly  enhanced  in  the  system  H2/CU.  Using  the  so-called  vibrational  efficacy, 


X  = 


A  Ev 

hUJy ib 


(32) 


where  A Ev  is  the  energetic  shift  between  the  sticking  curves  for  molecules  in  the  vibrationally 
ground  and  first-excited  state,  the  vibrationally  enhanced  dissociation  can  be  quantified.  The 
shift  A Ev  is  approximately  0.3  eV,  as  indicated  in  Fig.  8.  Consequently,  with  the  vibrational  fre¬ 
quency  of  H2,  fuvv ib  =  0.516  eV,  the  vibrational  efficacy  becomes  x  ~  0.6  which  is  often  interpreted 
to  mean  that  60%  of  the  vibrational  energy  is  used  to  overcome  the  barrier  for  dissociative  adsorp¬ 
tion. 


Vibrationally  enhanced  dissociation  has  been  known  for  years  in  gas  phase  dynamics  [71].  The 
basic  mechanism  can  be  discussed  within  a  two-dimensional  elbow  plot  as  shown  in  Fig.  6b.  The 
system  H2/Cu(100)  exhibits  a  barrier  which  is  located  after  the  curved  region  of  the  minimum 
energy  path  which  is  called  a  late  barrier.  For  such  a  topology  of  the  PES,  initial  vibrations  can 
be  very  helpful,  as  is  illustrated  with  two  typical  trajectories  in  Fig.  6b.  Initially  non-vibrating 
molecules  with  a  kinetic  energy  less  than  the  barrier  height  are  reflected  at  the  adsorption  barrier 
(dashed  line).  However,  if  the  molecule  is  already  initially  vibrating,  i.e.,  if  it  is  oscillating  back 
and  forth  in  the  d-direction,  then  the  vibrational  energy  can  be  very  efficiently  used  “to  make  it 
around  the  curve”  and  enter  the  dissociation  channel.  In  addition,  there  are  vibrationally  adiabatic 
effects  associated  with  the  lowering  of  the  vibrational  frequency  perpendicular  to  the  reaction  path 
which  also  contribute  to  the  vibrationally  enhanced  dissociation  [72]. 

In  addition  to  the  vibrational  state,  also  the  rotational  state  can  have  a  significant  influence  on  the 
dissociative  adsorption  probability.  Usually  additional  rotational  motion  suppresses  the  adsorption 
probability  because  the  molecule  will  rotate  out  of  a  favorable  orientation  for  dissociation  during 
the  adsorption  process  [8].  Invoking  the  principle  of  microscopic  reversibility  of  detailed  balance, 
this  leads  to  the  so-called  rotational  cooling  in  desorption,  i.e.,  the  mean  rotational  energy  is  below 
the  value  expected  for  thermal  equilibrium.  Conversely,  the  vibrationally  enhanced  dissociation 
corresponds  to  vibrational  heating  in  desorption. 


6  The  full  concert:  molecular  and  dissociative  adsorption  and 
scattering 


In  the  examples  discussed  so  far,  we  focused  either  on  the  energy  transfer  in  atomic  and  molecular 
adsorption  or  on  the  bond-breaking  process  in  dissociative  adsorption.  However,  there  are  molecule- 
surface  systems  in  which  all  these  processes  can  occur  at  the  same  time.  One  example  is  the 
interaction  of  02/Pt(lll).  This  is  a  well-studied  system  motivated  by  the  fact  that  the  adsorption 
of  oxygen  on  platinum  represents  one  of  the  fundamental  microscopic  reaction  steps  occuring  in  the 
car-exhaust  catalyst.  On  Pt(lll),  there  exist  both  physisorbed  and  chemisorbed  molecular  oxygen 
species  [73,  74]  as  well  as  dissociatively  adsorbed  oxygen.  The  chemisorbed  species  have  also  been 
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identified  in  total-energy  calculations  [75,  76]  using  density  functional  theory  (DFT)  within  the 
generalized  gradient  approximation  (GGA)  [77]. 

According  to  molecular  beam  experiments  [78,  79],  the  sticking  probability  of  02/Pt(lll)  first 
exhibits  a  strong  decrease,  and  then  after  passing  a  minimum  at  approximately  0.15  eV  the  sticking 
probability  levels  off  at  a  value  of  about  0.3  [78,  79].  Surprisingly,  in  molecular  beam  experiments  it 
was  also  found  that  oxygen  molecules  do  not  dissociate  at  cold  Pt  surfaces  below  100  K  [33,  80,  79], 
even  at  kinetic  energies  above  1  eV  which  are  much  greater  than  the  dissociation  barrier. 

The  theoretical  description  of  the  adsorption  dynamics  of  02/Pt(lll)  represents  a  significant 
challenge.  On  the  one  hand,  a  realistic  PES  is  needed  that  reliably  describes  both  the  molecular 
as  well  as  the  dissociative  adsorption  channels.  On  the  other  hand,  molecular  trapping  processes 
can  only  be  reproduced  if  the  energy  dissipation  to  the  platinum  substrate  is  properly  taken  into 
account,  as  discussed  in  sect.  4.  Direct  ab  initio  molecular  dynamics  simulations  represent  a  scheme 
that  meets  these  requirements,  but  it  is  computationally  still  very  expensive  [2,  81].  Using  empirical 
classical  potentials,  almost  arbitrarily  many  trajectories  can  be  computed,  however,  there  are  no 
reliable  interaction  potentials  available  treating  reactions  on  the  surface  and  the  surface  recoil  upon 
impact  on  an  equal  footing. 

As  an  intermediate  method,  tight-binding  molecular  dynamics  simulations  have  been  performed  [82, 
83]  with  the  parameters  of  the  tight-binding  Hamiltonian  derived  from  ab  initio  calculations  [84,  75, 
76].  In  tight-binding,  the  exact  many-body  Hamiltonian  is  replaced  by  parameterized  Hamiltonian 
matrix  elements  of  the  effective  one-particle  Hamiltonian  in  an  atomic-like  basis  set.  The  atomic- 
like  basis  functions  are  usually  not  considered  explicitly,  but  the  matrix  elements  are  assumed  to 
have  the  same  symmetry  properties  as  matrix  elements  between  atomic  states.  The  evaluation  of 
the  tight-binding  Hamiltonian  still  requires  the  diagonalization  of  a  matrix,  but  it  is  about  three 
orders  of  magnitude  faster  than  corresponding  ab  initio  calculations. 

A  comparison  between  the  calculated  [83]  and  the  measured  [78,  79]  sticking  probabilities  is  shown 
in  Fig.  9a.  It  is  obvious  that  the  experimental  data  are  qualitatively  and  even  semi-quantitatively 
reproduced  by  the  ab  initio  based  tight-binding  molecular  dynamics  calculations.  It  is  important 
to  note  that  also  the  experimental  finding  that  O2  does  not  directly  dissociate  upon  adsorption  is 
confirmed  by  the  calculations. 

There  is  a  simple  explanation  for  this  result  in  terms  of  the  topology  of  the  elbow  plots  (see  Fig.  9b). 
Dissociation  corresponds  to  an  event  in  which  the  molecules  enter  the  exit  channel  towards  the 
lower  right  corner  of  the  figures.  To  enter  this  channel  directly  from  the  gas  phase  through  the 
molecular  adsorption  state  requires  a  sharp  turn  of  the  trajectories.  In  Fig.  9b,  a  trajectory  of  an 
O2  molecule  directly  aimed  at  the  molecular  precursor  state  is  included.  Its  kinetic  energy  of  0.6  eV 
is  much  higher  than  the  dissociation  barrier  (~  0.2  eV  [76]);  still  it  does  not  dissociate.  It  becomes 
accelerated  by  the  attractive  potential,  hits  the  repulsive  wall  of  the  potential  and  is  scattered 
back.  This  shows  that  direct  dissociation  of  02/Pt(lll)  is  not  impossible,  but  it  is  very  unlikely. 
Thus  it  follows  that  because  of  this  steric  hindrance,  dissociation  of  O2  on  Pt(lll)  is  a  two-step 
process.  First  the  molecule  becomes  trapped  and  accommodated  in  the  molecular  chemisorption 
state,  and  only  subsequently  it  dissociates  at  sufficiently  high  surface  temperatures  due  to  thermal 
fluctuations  which  will  make  the  O2  molecules  enter  the  dissociation  channel. 
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Figure  9:  O2  adsorption  on  Pt(lll).  a)  Trapping  probability  of  02/Pt(lll)  as  a  function  of  the 
normal  component  of  the  kinetic  energy.  Results  of  molecular  beam  experiments  for  surface  tem¬ 
peratures  of  90  K  and  200  K  (Luntz  et  al.  [78])  and  77  K  (Nolan  et  al.  [79])  are  compared  to 
tight-binding  molecular  dynamics  simulations  for  the  surface  initially  at  rest  ( Ts  =  0  I().  b  and 
c)  Elbow  plots  of  the  PES  of  the  dissociation  of  02/Pt(lll)  determined  by  the  ab  initio  derived 
tight-binding  Hamiltonian  [76,  83].  The  configurations  of  the  remaining  O2  degrees  of  freedom  are 
illustrated  in  the  insets.  The  contour  spacing  is  0.2  eV  per  O2  molecule.  In  (b)  a  trajectory  of  an 
O2  molecule  with  an  initial  kinetic  energy  of  0.6  eV  scattered  at  Pt(lll)  is  also  plotted. 


The  molecular  dynamics  simulations  also  showed  that  contrary  to  common  belief  [33,  79]  the 
strong  initial  decrease  of  the  sticking  probability  is  not  caused  by  the  trapping  into  a  shallow 
physisorption  state.  Instead,  the  high  sticking  probability  at  low  kinetic  energies  is  again  caused 
by  the  steering  effect  which  becomes  quickly  suppressed  for  higher  kinetic  energies.  Thus  it  is 
not  the  energy  transfer  per  se  that  determines  the  sticking  probability  at  low  kinetic  energies  but 
rather  the  probability  to  enter  the  molecular  chemisorption  state.  All  molecules  that  find  their  way 
to  the  molecular  chemisorption  state  at  low  kinetic  energies  do  in  fact  remain  trapped. 

The  magnitude  of  the  steering  effect  is  a  consequence  of  the  strong  corrugation  and  anisotropy  of 
the  PES  of  02/Pt(lll).  The  significant  corrugation  is  illustrated  in  Figs.  9b  and  c.  The  lateral 
position  of  the  O2  center  of  mass  is  only  shifted  by  about  1 A  between  the  two  cuts.  Nevertheless, 
there  is  no  longer  any  chemisorption  well  present  but  rather  a  large  barrier  of  about  1  eV  towards 
dissociative  adsorption  which  becomes  even  larger  for  the  molecule  directly  at  the  on-top  site.  In 
fact,  the  majority  of  adsorption  pathways  are  hindered  by  barriers;  direct  non-activated  access  to 
the  adsorption  states  is  possible  for  only  a  small  fraction  of  initial  conditions. 

The  strong  corrugation  and  anisotropy  of  the  PES  has  further  consequences.  First  of  all  it  leads 
to  a  high  probability  of  dynamical  trapping  due  to  the  conversion  of  the  initial  kinetic  energy 
into  internal  molecular  degrees  of  freedom  which  for  this  particular  system  is  almost  independent 
from  the  kinetic  energy.  This  causes  the  leveling  off  of  the  trapping  probability  at  higher  kinetic 
energies.  Furthermore,  in  Fig  9a  the  results  of  calculations  for  non-normal  incidence  are  shown 
for  two  different  total  kinetic  energies,  Ei  =  0.3  eV  and  1.1  eV.  Additional  parallel  momentum 
strongly  suppresses  the  sticking  probability,  in  agreement  with  the  experiment  [78].  It  is  obvious 


RTO-EN-AVT-142 


4-19 


Simulation  of  Gas-Surface  Dynamical  Interactions 


ORGANIZATION 


Figure  10:  Angular  distribution  of  O2  scattered  from  Pt(lll)  with  an  angle  of  incidence  of  60°.  a) 
Calculated  angular  distribution  (surface  temperature  Ts  =  300  K)  for  an  initial  rotational  energy  of 
0.1  eV  (filled  circles)  and  0.3  eV  (open  circles),  b)  Comparison  of  the  measured  and  the  calculated 
angular  distribution  in  in-plane  scattering  for  an  angle  of  incidence  of  60° .  The  initial  kinetic  energy 
in  the  experiment  was  1.27  eV  [85]  while  the  TBMD  simulations  have  been  performed  for  an  initial 
energy  of  1.1  eV.  Theoretical  results  are  shown  for  both  an  initial  rotational  energy  of  0.1  eV  and 
0.3  eV. 

that  the  trapping  of  02/Pt(lll)  does  not  obey  normal  energy  scaling,  i.e. ,  it  is  not  a  function  of 
the  normal  kinetic  energy  alone,  but  total  energy  scaling  is  also  not  obeyed  since  the  results  for 
non-nornral  incidence  are  smaller  than  those  for  normal  incidence. 

For  non-nornral  incidence,  also  the  angular  distribution  of  the  scattered  O2  molecules  that  is  shown 
in  Fig.  10a  has  been  analyzed  in  detail.  The  distribution  shows  that  there  is  predominantly  in¬ 
plane  scattering,  i.e.  the  molecules  do  not  significantly  change  their  azimuthal  angle.  In  Fig.  10a, 
the  results  for  the  inplane  scattering  are  compared  to  the  experiment  [85] .  There  is  a  rather  good 
agreement  between  theory  and  experiment. 

However,  it  is  important  to  note  that  the  angular  distribution  of  Ar  atoms  scattered  from  Pt(lll) 
is  very  similar  to  that  of  O2  scattered  at  the  same  surface  [86,  85].  The  interaction  potential  of 
noble  gas  atoms  with  low-index  metal  surfaces  is  usually  relatively  structureless  and  only  weakly 
corrugated.  The  similarity  between  Ar/Pt(lll)  and  02/Pt(lll)  scattering  indicates  that  also 
the  O2  scattering  corresponds  to  the  reflection  from  a  rather  flat  surface.  This  has  already  been 
indicated  by  the  weak  out-of-plane  scattering  in  the  system  02/Pt(lll).  The  width  of  the  in¬ 
plane  distribution  can  in  fact  be  explained  by  the  energy  transfer  to  a  vibrating  flat  surface  [85] 
which  leads  to  a  certain  width  in  the  distribution  of  the  normal  component  of  the  kinetic  energy. 
Together  with  the  conservation  of  the  parallel  momentum  of  the  scattered  particles  this  causes  the 
broadened  angular  distribution  around  the  specular  direction  in  the  scattering. 

Thus  there  is  the  seemingly  paradox  conclusion  that  adsorption  experiments  suggest  that  the 
O2 /Pt (11 1)  interaction  potential  should  be  strongly  corrugated  while  scattering  experiments  in¬ 
dicate  a  rather  small  corrugation.  This  contradiction  is  caused  by  the  fact  that  the  scattered 
molecules  are  reflected  directly  from  the  repulsive  tails  of  the  potential  which  is  less  strongly  cor¬ 
rugated  while  in  the  trapping  the  corrugation  and  anisotropy  of  the  potential  energy  surface  closer 
to  the  surface  become  important. 
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Figure  11:  Illustration  of  the  different  apparent  corrugation  of  a  surface  as  a  function  of  the  angle 
of  incidence. 

This  phenomenon  is  illustrated  in  Fig.  11.  In  the  system  02/Pt(lll)  scattering  corresponds  to  the 
majority  channel  at  higher  kinetic  energies  while  adsorption  is  the  minority  channel.  Adsorption  is 
strongly  dependent  on  the  corrugation  of  the  surface,  however,  for  the  scattered  molecules  at  higher 
angles  of  incidence,  due  to  shadowing  effects  the  surface  looks  rather  smooth.  Since  adsorption  is 
the  minority  channel,  the  additional  scattering  flux  which  results  from  the  suppression  of  adsorption 
at  higher  angles  of  incidence  does  not  crucially  influence  the  scattering  distribution. 


7  Conclusions 


In  this  chapter,  I  have  briefly  reviewed  the  fundamentals  of  the  dynamics  of  gas-surface  interactions 
and  their  theoretical  treatment.  There  are  some  phenomena  in  this  interaction  such  as  elastic 
scattering  and  diffraction  that  can  only  be  understood  in  a  quantum  treatment.  Most  of  the 
processes,  however,  can  be  understood  qualitatively  and  even  quantitatively  based  on  classical 
dynamics  simulations. 

While  atomic  and  molecular  adsorption  require  a  appropriate  treatment  of  the  excitation  spectrum 
of  the  semi-infinite  solid  in  order  to  correctly  describe  the  energy  transfer  to  the  substrate,  in  the 
dissociative  adsorption  of  light  molecules,  in  particular  H2,  this  transfer  can  be  neglected  since 
the  crucial  process  is  the  bond-breaking  upon  adsorption.  However,  a  reliable  high-dimensional 
potential  energy  surface  is  needed  in  order  to  reproduce  the  dependence  of  the  adsorption  prob¬ 
ability  on  the  internal  degrees  of  freedom  of  the  impinging  molecules.  A b  initio  based  molecular 
dynamics  simulations  are  capable  of  giving  a  reliable  and  rather  complete  picture  of  the  interaction 
of  molecules  with  surfaces. 

Electronically  non-adiabatic  processes  have  not  been  discussed  in  this  chapter.  Their  theoretical 
treatment  is  still  problematic.  Such  processes  can  either  lead  to  further  energy  losses  due  to  the 
excitation  of  electron-hole  pairs  in  the  substrate,  or  they  modify  the  reaction  dynamics  due  to 
transitions  between  electronically  excited  states.  However,  for  most  thermal  processes  at  surfaces, 
electronically  non-adiabatic  processes  only  seem  to  play  a  minor  role. 
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